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Stability of a half-quantum vortex (HQV) in superfluid 3 He has been discussed recently by Kawakami, 
Tsutsumi and Machida in Phys. Rev. B 79, 092506 (2009). We further extend this work here and consider 
the A 2 phase of superfluid 3 He confined in thin slab geometry and analyze the HQV realized in this 
setting. Solutions of HQV and singly quantized singular vortex are evaluated numerically by solving the 
Ginzburg-Landau (GL) equation and respective first critical angular velocities are obtained by employing 
these solutions. We show that the HQV in the A 2 phase is stable near the boundary between the A 2 and 
Ai phases. It is found that temperature and magnetic field must be fixed first in the stable region and 
subsequently the angular velocity of the system should be increased from zero to a sufficiently large value 
to create a HQV with sufficiently large probability. A HQV does not form if the system starts with a fixed 
angular velocity and subsequently the temperature is lowered down to the A 2 phase. It is estimated that the 
external magnetic field with strength on the order of 1 T is required to have a sufficiently large domain in 
the temperature-magnetic field phase diagram to have a stable HQV. 
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1. Introduction 

A half-quantum vortex (HQV) in the A phase of superfluid 
3 He has been proposed first by Volovik and Mineev in 1976. 1 
They fully utilized the peculiar structure of the order param- 
eter manifold M = S 2 x SO(3)/Z 2 of the A phase, where 
S 2 is a manifold where the unit magnetic vector d resides, 
while S 0(3) is a group manifold representing the orbital de- 
grees of freedom fh+ih. They are intertwined by the "Z 2 factor, 
showing the uniqueness of the magnetic part and orbital part 
are required only up to the sign flip; pairs (d, m + in) and 
(-d, -(in + in)) define the same ordered state. 

Since then, there have been several theoretical 2 ^* as well as 
experimental 5 ' 6 works devoted to the stability of a HQV in a 
thin film of superfluid 3 He-A. Recently, Kawakami, Tsutsumi 
and Machida published papers studying the Majorana modes 
trapped in a HQV, 7 ' 8 inspired by the NMR experiment con- 
ducted by Yamashita et al. 9 They show that a HQV is stabi- 
lized when the system is rotated and obtain the phase diagram 
in rotation velocity-system size space. To date, however, the 
existence of a HQV has not been demonstrated in superfluid 
3 He yet in spite of intensive challenges. 

In the present paper, we further pursue their work and show 
that a stable HQV exists in the A2 phase near the boundary 
between the Ai and the A2 phases. The stability domain of a 
HQV has the temperature range comparable to that of the A 1 
phase if the exernal magnetic field has a strength on the order 
of 1 T. 

The rest of the paper is organized as follows. In Sec. II, we 
introduce the Ginzburg-Landau free energy of superfluid 3 He 
to establish our notation and convention. The order parameter 
describing a HQV is introduced in Sec. III. The critical an- 
gular velocity is defined in Sec. IV for a HQV and a singular 
vortex. Section V is devoted to numerical analysis. The order 
parameter profiles and comparison between the formation en- 



ergy of a HQV and a singular vortex are given. Section VI 
concludes this paper. 

2. Ginzburg-Landau Free Energy 

Let A a i be the order parameter of a superfluid phase of 3 He, 
where a is the spinor index while i is the orbital index. Then 
the bulk free energy in the Ginzburg-Landau expansion takes 
the form 

F B = - aA* ai A ai + PiA^fApjApj 

+ P2A* ai A ai A*fsjA pj + faAl^AvjAfsj 

+ faKA^Aj + feKi A l}i A f})Kp ( 1 ) 

where the coefficient a of the second order term has a temper- 
ature dependence a — a't with a constant a' and t = 1 - T/T c . 
We take account of the effect of the strong coupling through 
the paramagnon parameter 5 in the fourth order terms f3j as 

ft = -(1 + 0.l5)f3 , j3 2 = (2 + Q.28)p Q , 

y8 3 = (2 - 0.05<5)A), jS 4 = (2 - 0.555)jSo, 

#s = -(2 + 0.7<5)A). (2) 

The gradient free energy is given by 

F c =KidiA aj diA* aj + K 2 d i A ai d j A' aj 

+ K 3 diA aj djA: yi . (3) 

The coefficients K/ satisfy 

K 1 =K 2 = K 3 =K (4) 

in the weak coupling limit. We emply the relation (4) in the 
rest of this paper for simplicity. As a result, the coherence 
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length is uniquely defined as 




(5) 



'K 1 

In the following, we take superfluid 3 He confined between 
two parallel plates, the distance of which is less than the 
dipole coherence length. The superfluid is rotated around the 
z-axis, which is perpendicular to the plates, and subject to a 
strong magnetic field along the z-axis so that the superfluid is 
in the A2 phase. The parallel plates introduces the boundary 
condition such that the / vector is perpendicular to the plates 
at the boundary. As a result, the orbital states of the order pa- 
rameter are restricted to l z = 1 or L = - 1 at the boundary. A 
strong magnetic field along the z-axis aligns the d vector in 
the xy-plane. 

It turns out to be convenient in this setting to change the 
basis of the order parameter from (x,y,z) to (1,0,-1), char- 
acterizing the z-component of the angular momentum. These 
two sets of basis vectors are related as 



V2 V 



± w 



eo 



e 7 . 



(6) 



From now on, we change the notation of the order parameter 
from A ai with respect to the e, = (e x , e y , e z ) basis to A^ v with 
respect to e v = (c_,Co,c+) basis. 

The magnetic field coupled to the superfluid changes the 
second order term in F B to 

- £ [(a + i!H)A* +v A +v + (a- r}H)A*-A-y] 

v=± 

= - J] a't[(l + h)A' +v A +v + (1 - ftMlyA-y], (7) 



where h = h/t. The variable h = r]H/a' is a dimensionless 
parameter corresponding to the magnetic field strength. The 
parameter r/ is a constant yielding coupling between H and 
the condensate. 

Let us analyze a uniform superfluind in the A2 phase with 
L = 1 by employing these free energies. Since A + _ = A__ = 
for this state, the bulk free energy reduces to 

F B =-a't(l+h)\A ++ \ 2 +{3 24 \A ++ \ 4 

-a'f(l-/0|A_ + | 2 +/? 24 |A_ + | 4 

+ 206 2 4 + 2 y 6 5 )|A ++ | 2 |A_ + | 2 , (8) 

where p 24 = p 2 + Pa = (4 - 0355)p , P24 + 2/3 5 = -1.75(5$). 
In the weak coupling limit 6 = 0, we obtain p 24 + - 
and as a result A ++ and A_ + decouple. We also note that in 
case 6 > 0, we obtain /?24 + 2/?5 < and hence the coupling 
between A ++ and A^+ is attractive. 

We introduce the following rescaling of physical quantities 
to simplify the notations; 

Pi -» pffii, |A ++ | 2 , |A_ + | 2 -> ^|A ++ | 2 , ^|A_ + | 2 , 

Po Po 



F -, (a ' t) F 
r B — > — r B - 

Po 



(9) 



Then the bulk free energy becomes 

F B = -{l+h)\A ++ \ 2 +p 2A \A, 



-(l-/0|A_ + | 2 +p- 2 4|A- + | 4 
+ 206 2 4 + 2^5)|A ++ | 2 |A_ + | 2 . 



(10) 



The bulk order parameter is fixed by minimizing Eq. ( 1 0) with 
respect to |A ++ | 2 and |A_+| 2 as 



A ++ | 2 



(1 +h)p 2A -(\ -h)(J3 2 4 + 2p 5 ) 



\A- + 



-8p 5 (J3 2A + p 5 ) 

2 _ (l-h)p 24 -(l+h)(p24 



2p 5 ) 



(11) 



-&p 5 (p 24 +p 5 ) 

We take, without loss of generality, the following sign con 
vention 



l(0) 



A_ + =A<°> 



(l+h)p 2 ,-(l-h)(J3 2 , + 2p 5 ) 
\ " -8p 5 (J3 24 + p 5 ) 

l(l-h)p 24 -(l+h)(p 24 +2p 5 ) 



-8p 5 (p 24 +p 5 ) 



(12) 



This choice gives the {/-vector parallel to the x-axis in the A 
phase resulting in the limit h — > 0. 

The coherence lengths of A ++ and A_ + in the presence of 

h * are ^ + = 1/ +h) and = 1/^(1 -h) and they 
satisfy inequalities 



? + < ? < t- 



(13) 



We now look at the gradient energy Fq- Consider a vortex 
along the z-axis and assume the order parameter is transla- 
tionally invariant along this axis. Let us introduce the cylin- 
drical coordinates (r, <p,z) assuming d/dz is a null operator. 
Let n^y e Z be the quantum number of the component A^ 
and write it as 



(14) 



The gradient term introduces the coupling between the orbital 
components + and -, namely the coupling between A M+ and 
A^-. The quantum numbers n^ ± must satisfy the condition 



+ 2 



(15) 



for the vortex to be cylidrically symmetric around the z-axis. 
If this is the case, the gradient energy takes the form 



F G = J] 



—r - HOLr 
or r 



V —C^y 

or r 



■c 



IJ.V 



2 



dr " v 



(16) 



where the rescalings r — > %r and F G — > [(a'f) 2 /jSo] Fq have 
been made as before. 

3. Half-Quantum Vortex 

The order parameter of a HQV proposed by Volovik and 
Mineev 1 takes the form 

A ai = A A d a (m + ih)i 

= A A e i(fil (cos |c. t + sin |« y ) (e x + ie y )i (17) 
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in the A-phase with vanishing magentic field H — 0, where it 
is assumed that the /-vector is directed along the z-axis, while 
the rf-vector points in the xy-plane. Equation (17) is rewritten 

as 



A A (e + - e^e-)^ 



(18) 



This shows that the order parameter of the HQV represented 
in the (1,0,-1) basis has a non-vanishing winding number 
only in the component c_. Similarly, there is an order param- 
eter of a HQV, in which only the component e+ has a non- 
vanishing winding number. 

By considering the condition (15), the order parameter (18) 
yields a vortex with quantum numbers 



((n ++ , «+_), («-+, «_)) = ((0, 2), (1, 3)) . 



(19) 



We call this vortex as a vortex of type (0, 1) to distinguish 
it from other types of vortices introduced in the following. 
When the superfluid is rotated in the opposite sense, the re- 
sulting vortex has an order parameter in which e ,ip is replaced 
by e _!(0 in Eq. (18), which will be called a vortex of type 
(0,-1) having quantum numbers 



((« ++ , »+_), (n_ + , n_)) = ((0, 2), (-1,1)). 



(20) 



It is important to realize that the structure of a vortex of type 
(0, 1), obtained by rotating the superfluid in the positive sense 
with respect to the /-vector, is different from that of a vortex 
of type (0,-1) obtained by rotating the superfluid in the opp- 
site direction. The condensate with orbital angular momentum 
spontaneously breaks the rotational invariance and hence the 
clockwise rotation and anticlockwise rotation are not mirror 
reflections of each other. 

4. First Critical Angular Velocity 

Let R be the radius of a cylindrical container and Q be the 
angular velocity with which the cylinder rotates. Now we ob- 
tain the condition under which a vortex stably exists at the 
center of the container. The gradient free energy in the rotat- 
ing system is obtained by replacing the (^-derivative as 



dm cL 2m 

— -> — - -r-(n X #%, 

ir ir n 



(21) 



where m is the mass of a 3 He atom. Let us first consider a 
HQV, in which the component e_ with lower creation en- 
ergy has a non-vanishing quantum number n_+. There are two 
terms of the form 

1 2m \ 2 1 

-Qr C_+(r) 2 (22) 

r n J 

in the gradient energy (16). The coefficient of a term linear in 
Q is nothing but the angular momentum and the total angular 
momentum of the system is found to be 

Z,H = 2 x An— [ R rdrC- + (r) 2 = An— (A (0 |) 2 R 1 , (23) 

% Jo h v ' 

where we noted that the contribution of the vortex core to the 
total angular momentum is negligible. 

The vortex formation energy measured with respect to the 
uniform bulk energy Fo is evaluated as 

Fvor = 2tt J rMF - Fo) = An (a™ f QnR + C-). (24) 



The parameter C will be evaluated numerically later. The 
first critical angular velocity for a formation of a vortex with 
spin component, namely a vortex in the spin component H 
is obtained by solving 



as 



n 



r(-) 



QL = 



n 



2mR 2 



QnR + C) 



(25) 



(26) 



2m( A (oyf R2 

An angular velocity will be scaled by h/2mR 2 from now on. 
As a result, the crtical angular velocity is written as 

Q< _) = InR + C-. (27) 

A singular votex (SV) with a winding number 1 is obtained 
by setting the quantum numbers of A MV to 

((n ++ ,n + _),(n_ + ,n__)) = ((1,3), (1,3)). (28) 

The total angular momentum of a S V is 

L (s) =2 x 4?r^p ^ rdr [c ++ {rf + C_ + (r) 2 ] 



R 1 



The formation energy of a SV is 



■ An 



(A<°>) 2 + (A?») 2 



(In/? + C,) 



and the first critical angular velocity is 
Qi s) = lnfl + C s 



(29) 



(30) 



(31) 



in the dimensionless form. 

Whether a HQV forms or a SV forms as the angular veloc- 
ity increases depends on the magnitude relation between 
and 0* ?) . Let us consider the case 6-0 and h > to begin 
with. It follows from the inequality £ + < that a vortex in 
A_ + is energetically favorable than that in A ++ and it is ex- 
pected that _) < Q* is satisfied. In case 5 > and h = 0, 
the coupling between A ++ and A_ + is attractive and a SV is 
expected to be favorable compared to a HQV. This is because 
low magnitude A ++ and A_ + overlap at the common vortex 
core in a SV while they do not in a HQV, thus gaining more 
negative energy for the former. Then an inequality Cl^ > 
is expected to be satisfied. 

It is expected from the above arguments that a S V is formed 
first as the angular velocity is raised from zero when 6 > and 
h is small. When the external magnetic field is strong enough, 
in contrast, there is a region in the temperature-angular veloc- 
ity domain in which a HQV is formed first. These statements 
will be verified numerically in the next section. 

5. Numerical Analysis 

We have solved the Ginzburg-Landau equation with respect 
to (C++(r), C+-{r), C-+(r), C—{r)) numerically. Four choices 
of the quantum numbers are considered; 

(a) A HQV (0,1) with ((n ++ , n + _), (n_+, «__)) 
((0,2), (1,3)), 

(b) A SV (1) with (( B++ , n + _), («-+, n__)) = ((1,3), (1,3)), 
both with Q. > and 
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-1) with ((«++,«+_),(«_+ 

with ((«++, «+-), («-+, « )) 



(c) A HQV (0,-1) with ((n ++ ,n + _), (n_ + ,n__)) 
((0,2), (-1,1)), 

(d) A SV (-1) 
((-1,1), (-1,1)), 

both with Q < 0. The parameters h and 5 are changed from 
0.0 to 0.5 with a step 0.1. The boundary condition at r — R 
does not affect the formation energy since R » 1 is assumed. 
We take the boundary condition 



(A ++ (R),A+_(R),A_ + (R),A__(R)) 



(32) 



which corresponds to a vortex embedded in a uniform / tex- 
ture with L = + 1 . 

The order parameter profiles for HQV (0, 1), SV (1), 
HQV (0,-1) and SV (-1) with 5 = 0.2, h = 0.2, and R = 30 
are shown in Fig. 1 . 

Next the first critical angular velocities are obtained by 
evaluating the free energies of HQV's and SV's with our nu- 
merical solutions and then employing Eqs. (27) and (31). Fur- 
thermore, we repeat the same calculation with R — 40 and 50 
and fit the first critical angular velocities thus obtained with a 
function £l c — A InR + C, A and C being constants. The result 
shows that A is in fact 1 with a good precision, as expected, 
and we have determined the h- and ^-dependences of C and 
C s , the C value of a HQV (0, 1) and a SV (1), respectively. 
Figure 2 shows C - Q. c -InR and C s = Q. c -InR for cases 
(a) and (b) as functions of h for 5 = 0.2. The dimensionless 
magnetic field h also takes values h = 0.0,0.1, . . .,0.5. The 
result shows that, for 5 + 0, there exists h c at which the in- 
equality > Q* s) flips to £^ _) < Q. { c ] as h is increased. The 
critical magnetic field h c vanishes for 6 = 0, showing there is 
a range of Q in which a HQV is stable for any h > 0. The first 
critical angular velocity h c has been estimated by finding the 
intersection of numerically interpolated curves O^T and Oj. !) 
as functions of h. Figure 3 depicts the ^-dependence of h c thus 
obtained. In case we reach a region in the Q. axis 

in which a SV is stablized if the angular velocity is further 
increased beyond £l c . The boundary between two stability 
regions along the Q-axis is found from 



(33) 



(-)->« 



- In R. The region in 



Figure 2 also shows C(-)_>(» = O!, 
which a HQV is stable is bounded by two curves an d 
Ojr ) in Fig. 2. 

Next, the /i-dependences of C and C s for cases (c) and (d), 
respectively, are depicted in Fig. 4. They are different from 
those of cases (a) and (c), reflecting upon the difference in the 
vortex structures for Q > and Q < 0. Figure 4 also shows 
fl[ ' , similarly to Fig. 2. A HQV is stable in the region 
bounded by two curves ~'~*( s > and £2^ ~ . Let us evaluate the 
critical value h c , at which the HQV-stable region appears. The 
critical value will turn out to be almost the same as that with 
Q > 0, in contrast with C and C s . The region in the ftf-plane 
(the temperature-pressure plane) where a HQV is stable is ob- 
tained from the ^-dependence of h c . The phase diagram for a 
given h = rjH/a' is shown in Fig. 5, where t — h/h has been 
used. The Ai phase is also shown in Fig. 5 for comparison. 
Although we have ignored the ^-dependences of rj and the su- 
perfluid transition temperature T c at H = 0, the comparison 
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Fig. 1. (Color online) Vortex solutions of the Ginzburg-Landau equation. 
Order parameters C++(r) are shown for r5 = 0.2, h = 0.2 and R = 30. Quan- 
tum numbers ((«++, «+-), («-+,« )) are (a) ((0, 2), (1, 3)), corresponding 

to a HQV with CI > (b) ((1,3), (1,3)), corresponding to a SV with 
f2 > (c) (0,2), (-1, 1)), corresponding to a HQV with Q < and (d) 
((-1, 1), (-1, 1)) corresponding to a SV with SI < 0. 



between the width of the Ai phase and that of the HQV-stable 
region is meaningful for a fixed 5. Figure 5 shows that the 
width of the HQV-stable region along the f//i-axis is compa- 
rable to that of the Aj phase. Morevoer the former increases 
compared with the latter for small 6 (low pressure) region. 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 5 



o - 

0.1 



-0.1 
-0.2 



InR 



C s 



sv 



V 



0.2 



(-)-►(») 



0.3 0.4 0.5 
HQV 



Meissner State 
8=0.2 



Fig. 2. (Color online) HQV-stable region for fl > 0. The solid line is the 
boundary between the Meissner state without a vortex and a HQV while 
the dotted line is the boundary between the Meissner state and a SV. For 
sufficiently large CI, transition from a SV to a HQV takes place as h is 
increased and the dashed line is crossed. The area bounded by the solid 
line and the dashed line is the region where a HQV has the least energy. A 
HQV-stable angular velocity region exists only for h > h c {5) when S > 0, 
while it exists for any h when 6 = 0. The intersecting point of the solid, the 
dashed and the dotted lines gives h c (S). 




Fig. 3. (Color online) HQV-stable region in the h-S domain for £1 > 0. The 
boundary corresponds to h c (S). Here HQV denotes the parameter region in 
the h-S plane where there exist £1 > for which a HQV has lower energy 
than a SV. There is essentially no difference in the graph for £1 < in the 
h-S plane. 
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Fig. 4. (Color online) HQV-stable region for fi < 0. The solid line, the 
dashed line and the dotted line denote the same boundaries as in Fig. 2. 



6. Conclusion and Discussion 

We have obtained conditions with which a half-quantum 
vortex stably exists and have shown that the stability region 
of a HQV has a comparable range to that of the A; phase 
along the t/h-axis. 

To obtain the HQV, the temperature and the magnetic field 




Fig. 5. (Color online) HQV-stable region in the t-S plane. The HQV-stable 
region is determined by the 5-dependence of h c . For comparison, the Ai 
phase, given h = riH/a', is also shown here. The origin of the horizontal 
axis corresponds to T = T c at H = 0. The critical temperature T c of the A[ - 
phase is 1 in the present scaling. The 5-dependence of r\ and T c at H = 
is ignored and, hence, comparison of the phases with different S should 
not be taken seriously. In spite of this, comparison of the widths of the A[ 
phase and the HQV-stable region is meaningful. It shows that the width 
of the HQV-stable region has the width of the same order of that of the 
A i phase. The width of the former increases as S is lowered (low pressure 
region). 



must be fixed in the HQV-stable region with no rotation first 
and subsequently the angular velocity must be increased be- 
yond the critical angular velocity. The opposite scenario, in 
which the system is rotated beyond the critical angular ve- 
locity first and then the temperature is lowered to form the 
A2 phase through the Ai phase, does not lead to a HQV 
formation. This is because a sigular vortex forms in the Tf- 
component while the system is in the Ai phase and it is impos- 
sible to eliminate this singular vortex after the system reaches 
the A2 phase. We believe that a magnetic field on the order 
of 1 T is required to have a large enough HQV-stable region 
comparable to that of the Ai phase. 

It is desirable to observe the direct NMR signal from a 
HQV for its detection. Nonetheless, direct observation can 
be rather challenging in the presence of a strong magnetic 
field. Note, however, that, when HQV formation takes place, 
there are two transitions associated with vortex formation in 
the vicinities of fi' an d as the angular velocity is in- 

creased. We expect these two transtions can be experimentally 
observable by one way or another. 
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